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SECTION  I 
INTRODUCTION 

With  reference  to  Fig.  I,  consider  the  flow  about  a finite-span  pure^- 

2 

jet-flapped  wing  with  partial  or  full-span  blowing  of  the  jet  sheets. 

The  motivation  for  study  of  this  flow  is  the  possible  use  of  jet  flaps 
for  maneuvering  combat  aircraft  at  high  subsonic  and  transonic  speeds. 

This  contrasts  with  the  more  conventional  application  of  the  jet-flap 
supercirculation  principle  wherein  the  jet  sheets  are  employed  to  augment 
the  lift  of  a mechanical  flap  during  the  takeoff  or  landing  flight  phase 
of  an  aircraft.  In  this  latter  application,  compressibility  effects  are 
of  secondary  importance. 

On  the  assumption  of  small  flow  perturbations  and  a restriction  to 
wings  without  dihedral  and  twist,  linear-subsonic  and  nonlinear-transonic 
similarity  rules  are  derived  herein  for  the  subject  flow.  In  developing 
these  similarity  rules  a new  jet-sheet  compatibility  condition,  which  is 
second  order  relative  to  the  jet-sheet  internal  flow,  is  derived.  The 
new  compatibility  condition  yields  a jet-momentum  coefficient  similarity 
parameter  differing  from  the  conventional  parameter  in  that  it  includes 
the  effect  of  the  jet-supply  pressure  ratio.  In  application  of  the 
similarity  laws,  particular  attention  is  given  to  camber-line  effects 
which  heretofore  have  been  unimportant  in  applications  involving  jet- 
augmented  mechanical  flaps. 


The  nonlinear  transonic  similarity  rules  derived  herein  for  jet- 
flapped  wings  are  new  (to  the  best  of  the  author's  knowledge).  For 
the  linear-subsonic  case,  however,  similarity  laws  for  jet-flapped  wings 
previously  have  been  presented  by  Siestrunk  (Ref.  1),  Levinsky  (Ref.  2), 
and  Elzweig  (Ref.  3).  Each  of  these  investigators  employs  a different 


^That  is,  a wing  employing  jet  sheets  alone,  unassisted  by  mechanical 
devices  such  as  flaps  or  ailerons. 

2 

For  convenience  in  discussion,  reference  to  a "blown"  or  “unblown"  wing 
will  connote,  respectively,  a wing  with  or  without  the  jet  flaps 
operating. 
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scaling  law,  and  only  Levinsky  considers  the  finite-span  wing;  the  other 
treatments  are  for  two-dimensional  flow.  The  present  treatment  is  more 
general  than  the  foregoing  ones  in  that  it  applies  to  both  two-  and  three- 
dimensional  wings,  allows  for  the  selection  of  the  scaling  law  most 
appropriate  to  the  particular  problem  under  consideration,  and  includes 
a second-order  scaling  of  the  jet  internal  flow  that  accounts  for  the 
jet-supply  pressure  ratio.  The  laws  employed  herein  reduce  to  those  of 
the  aforementioned  investigators  upon  appropriate  selection  of  the 
parameters. 

As  a consequence  of  the  assumption  of  small  perturbations,  application 
of  the  present  results  obviously  is  limited  to  small  thickness  and 
camber  ratios,  small  angles  of  attack,  and  small  jet-deflection  angles. 


2 
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SECTION  II 
ANALYSIS 

For  an  untwisted  wing  with  camber  and  thickness,  but  without 
dihedral,  the  geometry  may  be  given  by 


<t) 


Si 


(tc)  a 


cr  a 


) 


z$  = oF  (x*,  y*;  gi 
where 

F = ± 2r  fT(x*,  y*;  gi(T))+  r1  y*‘-  f x* 


Ml) 


with  a specified  spanwise  blowing  distribution  given  by 


c-j  “ cj  (F*j> 

6,  = 64  (y*4> 


Jj  - -J  ^ j' 


} 


(2) 


The  g^  quantities  in  Eq.  1 are  the  nondimensional  geometric  parameters 
required  to  define  the  specific  wing  geometry  being  considered,  with 


Si 

Si 


(T)  (T) 

_ (t) 

R<*> 

"Sj  > »2  ’ 

(<)  (<)  M ^ . 

Srix  ’ 
„ (<> 

=*1 

*> 

CM 

g\  , 

1*1,  2)  3* 

(3) 


For  some  configurations,  some  of  the  g.^T\  9-j^  parameters  may  be 
identical  (e.g.,  g^T^=  g1  The  parameter  a may  be  selected  to  be 

any  one  of  the  parameters  Tp,  <r,  a,  0^  or  6^  depending  upon  the  type 
of  problem  under  consideration. 


For  definiteness,  it  is  of  interest  to  apply  the  foregoing  relations 
to  a specific  class  of  wings  designated  as  "trapezoidal."  A "trapezoidal 
wing"  is  defined  herein  as  a wing  (without  dihedral  or  twist)  for  which 
the  half-span  planform  is  a trapezoid  with  the  root  and  tip  chords 
parallel,  the  wing-section  geometries  affine  to  each  other  at  varying 
spanwise  stations,  and  the  section  characteristic  dimensions  (chord, 
thickness,  camber  height,  etc.)  vary  linearly  in  the  spanwise  direction. 
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For  this  class  of  wings,  the  functions  in  Eq.  1 have  the  form 


fT  = fT(x*,  y*;  A tan  A,  \Q,  Xt) 


f<  = fK(x*,  y*;  A tan  A,  Ac , A^) 


where 


,(t> 


) 


U) 


A tan  A 
A 


~:k 

gi  = = o;  i > 


g2  C 

Xt,  8W  - X 


>3 


Explicit  forms  of  the  and  f^  - functions  are  not  displayed 
since  they  serve  no  useful  purpose  here.  They  are  easily  derived, 
however,  if  required. 


Deflected  ailerons  and  flaps  may  be  treated  by  introducing  additional 
"camber-like"  functions. 

The  governing  partial  differential  equation  for  the  perturbation 
velocity  potential  (Ref.  5)  is 


where  for  linear  subsonic  compressible  flow  j = 0,  and  for  nonlinear 
transonic  flow  j = 1 . 

The  shock-wave  compatibility  condition  (Ref.  5)  is 

62  (A$x)2  + (Hy)2  + (Aij>z) 2 = k (is)(*Xl  + 4>x2)(A<Px)2  (7) 

where  the  subscripts  x,  y,  and  z denote  partial  derivatives,  and  A 
denotes  a jump  in  the  modified  quantity  such  that  if  "1"  and  "2"  denote, 
respectively,  upstream  and  downstream  conditions  relative  to  the  shock 

4*x = % - V etc- 
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If  the  lateral-flow  velocity  perturbations  are  assumed  to  be  neglible 
compared  to  the  vertical  and  longitudinal  ones,  the  jet-sheet  compatibility 
condition  (see  Appendix  A)  is 

(8) 


A (At)  = £ (A0jS)  K c-i 
U.  W Q 4 


where 


K*  = 1 


K = 1 - ^ (1- 
* 2y  p/p 


for 

for 


(P/P.  1 p/P*> 
(P/p . < P/Poo  < 5 


(9) 

(10) 


For  y = 1.4,  P/p*  = 1.893 


We  now  introduce  the  transformations 

X = Sxx,  y = syy,  z = szz,  <£  = 

U = s..U  , 8 = sR6  , (U  k)  = sk  (U  k) 

oo  U oo  ’ oo  P oo  ’ CO  ' K 00 


} 


04) 


Flow  similarity  is  achieved  by  determining  appropriate  values  of  the 
s-stretching  factors  as  dictated  by  the  governing  equations  and  boundary 


The  boundary  conditions  at  x = - 00  and  y = ± °°  are 

. 

<J>  = * <t>  =0 

x y z 

(ID 

The  wing  surface  boundary  condition  is 

(M)  = u a f (x*  v*-  ^ 

^3z;  _ a 3x*  ' 'x  * y * a,  a,  a; 

Z=U 

(12) 

The  jet-sheet  boundary  condition  is 

1 

Ojs  (x*,  y*,  0)  - 04>/3z)z=()/uoo 

(13) 

a 
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We  consider  first  the  constraints  imposed  on  the  stretching  factors 
by  the  governing  partial  differential  equation,  Eq.  6,  and  the  shock 
compatibility  condition,  Eq.  7.  Applying  the  transformations  of  Eq.  14 
yields  the  governing  differential  equation: 


sf  n 2 if*  + si  iii  . s£  if*  = , c ii  iff 

sfs . 00  3X7  s$  3y4  s*  32*  sks* 


and  for  the  shock  compatibility  condition 

SX  -2  - , SV  - 2 sf  ~ 5 

TzTz  e.  (a*x)2  + -4  (a*  > + -f  (a*  > 


i 16) 


Comparison  of  Eqs.  6 and  15  for  the  linear  and  nonlinear  problems  and 
Eqs.  7 and  16  for  the  nonlinear  problems  yields  the  following  conditions 
for  flow  similarity 

Linear 


sx  „ „2  „ _2 

-4  sy  sz 


s s s 

— ( ) = an  arbitrary  constant 

Sk  % 


Nonlinear 


»y 

o2  y 


For  the  linear  case,  the  freedom  in  choice  of  the  parametric 
3 

combination  ( susx  /s^)  introduces  the  option  of  an  additional  degree 
of  freedom  in  the  scaling.  This  option  is  a well-known  property  of 
subsonic  scaling  laws  (see,  e.g..  Refs.  5 and  6). 


Equations  17a  and  17b  may  be  generalized  in  the  following  manner 


S r8u!S 

e/e  ( sa  ’ 
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where 


e/e  » an  arbitrary  constant 

e/e  » (Dj/Uk) 

00  00 


(Linear)  (19) 

(Non-linear)  (20) 


For  the  linear  case,  the  ratio  (e/e)  is  analogous  to  the  ratio  (X/X) 

.3 


employed  by  Spreiter  (Ref.  5).  The  two  ratios  are  related  by 

(e/e)  = (e«,2x/e2  x). 


Equating  the  first  and  third  terms  in  Eq.  18  yields  s = sfls,,  from 

X p z 

which  it  may  be  determined  that 


(susX) 


/SuSzN 

se  5 


(21) 


In  the  analyses  which  follow,  (s^s^)  always  appears  in  one  of  the 


parenthetical  combinations  shown  in  Eq.  21.  The  stretching  factors  s 


and  s^  therefore  need  not  be  determined  separately  and  the  quantities 


(susx/s^  and  (s^/s^)  may  be  treated  as  unknown  parameters. 


Substituting  Eq.  21  in  Eq.  18  yields,  with  ft/fl  = s^/s^ 


sx 

X o 


2 

sx  sg 
(g/e)(G/fl) 


(22) 


The  relation  Eq.  22  defines  three  equations  relating  the  stretching 
factors.  This  set  of  equations  may  have  a variety  of  forms,  depending 
upon  the  algebraic  manipulations  performed.  The  forms  used  herein  are 


By  (efi/e sx 

(23) 

8Z  (eO/efi)1^*  sx 

(24) 

o/o  - (6» /e»)3(e/e) 

(25) 

The  symbol  X is  Spreiter's  X,  not  the  X used  herein 

7 
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where 


/'*  ~ (t)  ~ (*)  Tr  Kr  “ \ 

2*  = cfF  (x*,  y*;  gi  , gi  g-,  j ) 


F = ± l£  f}(X*,  y*;  gi(T))  + I1  fK  (x*.  y*;  gi(K>)-  fx*  {31) 


The  principal  linear  dimensions  of  the  original  and  transformed 
wings  are  related  by 

cr  = sxcr;  b = syb;  tr  = sztr 


yi  = 8vyT;  vo  = syyG 


yielding 


x*  = x*;  y*  = y*;  (eR)^3z*  = (efi)^^ z* 

. ..  * (27) 

yI  yI’  yo  yo 

The  transformation  relations,  Eq.  14,  and  the  similarity  conditions  of 
Eqs.  17a  and  17b  also  yield 

(tan  A)/3OT  = (tan  A)/BM  ^g) 

SooA  — 3o<>A  (29) 

Any  one  of  the  parameters  s , s or  s in  Eq.  26  may  be  assigned  an 

x y z 

arbitrary  value.  For  unit  values  of  s , s , and  s respectively,  Eq.  26 

x y z 

indicates  that  wings  of  identical  root  chord,  span,  or  root  thickness 
respectively  are  being  compared  in  the  original  and  transformed  spaces. 

A common  selection  is  sx  = 1-0 

To  achieve  manageable  results,  it  is  necessary  to  assume  that 
?T  = f , fK  = f^,  and  F = F.  That  is,  the  analysis  is  confined  to 
families  of  wings  defined  by  these  constraints.  The  wing  geometry  in  the 
transformed  space  then  becomes 


cj  * £j  (yj *) 

fij  = (yj*) 
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The  equality  of  the  f's  requires  that 

- (t)  (t)  - (ic)  (k)  ( 33) 

8i  “ 8i  i 81  = gi 

where  for  trapezoidal  wings  the  gi 's  are  defined  by  Eq  5.  The  equality 
of  the  F's,  in  combination  with  the  equality  of  the  f's,  requires  that 

Tr/a  - tr/a  (34) 

<r/o  = Kr/o  (35) 

a/5  = a/a  (35) 

For  trapezodial  wings,  the  invariances  specified  by  Eqs.  33  require 

that 

A tan  A = A tan  A . _ . 

(37) 

A - A , A-  A A - A (38) 

C c’  t t*  IC  IC  v °' 

It  is  apparent  that  the  transformations  of  Eq.  14  are  consistent  with  the 
conditions  of  Eqs.  37  and  38  which  are  a consequence  of  restricting  the 
analysis  to  a family  of  wings.  Also  note  that  Eq.  37  is  equivalent  to  the 
pair  of  Eqs.  28  and  29. 

The  boundary  condition  on  the  wing  in  the  transformed  space  is  given 


^^z-o  " ^oo  0 ax*  F(x**  > 


Applying  the  transformations  of  Eq.  14  to  the  boundary  condition  in  the 
original  space  yields 

"4  ^ ° k*  y*:  — ' l4<1 


Equations  39  and  40  are  similar  if  (o/o)  * s^/s  s ).  This  condition, 

$ u z 

an  invariant  of  the  transformation,  may  be  written  as  (ll/o)  * (il/o).  In 
applying  the  similarity  laws,  the  aforementioned  invariant  is  frequently 
employed  in  the  form 

(0/5)w  - <fi/o)w  (41) 
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where  the  value  of  m is  arbitrarily  selected  to  achieve  a particular 
desired  form  for  the  invariant  parameters  in  the  final  formation^  the 
similarity  rules.  Common  choices  are  (-1)  and  (2/3)  respectively  for 
subsonic  and  transonic  flows. 


Since  the  jet  sheet  is  a stream  surface,  Eq.  44  applies  to  the  jet-sheet 
slopes  0js  and  0js.  A basic  parameter  for  the  jet  sheet,  however,  is  the 
slope,  0.,  at  the  wing  trailing  edge.  Applying  Eq.  44  to  0.  and  0. 

J J J 

yields 

(n/Oj)40  - (ft/0j>(‘)  (45) 

where 

0j  “ - («  + 6j),  0j  * - (a  + 8j)  (46) 

Equation  46  may  be  written  as 


0,  - - 6,  (|-  + l),  0 - - 61  (f-  + 1) 

J j J 

For  o * 6j,  Eq.  36  yields 
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When  o is  selected  to  be  6.,  an  optional  form  for  Eq.  45  is.  therefore, 

J 

(fl/6)w  = (ft/ 6j)“  (o  - 6.,)  (4£ 

The  compatibility  condition  for  the  jet  sheet  in  the  original  space 
is  given  by  Eq.  8;  in  the  transformed  space  it  is 

, aT  - 30. 

= c c (49 

U„  3x  z*0  4 ( 3J  )K7rcj  V 

Applying  the  transformations  of  Eqs.  14  and  44  to  Eq.  8 yields 

(1/V  c 1 1 _ 

(l/sujff^d/sx)  V9x  z*0  ■ 4 sx  (s^/susz)  s*  hi  > KitCJ 


i /-Ml  _ C SZ  /°  j S \ ^ 

U<c  z=0  4 sx  3x  7fC^ 

The  ratio  (sz/sx)  may  be  obtained  from  Eq.  24.  Substituting  in 
Eq.  50  yields 

f #i.o  ■ f #>1/3<^)  ^ 


From  Eqs.  49  and  51 


(efi)1^3K7rcj  = (efi)1/3^ 


where,  from  Eqs.  25,  41,  45,  and  48,  the  ratio  (S/ft)  required  for 
Eq.  52  is  given  by  any  one  of  the  following  equations 


tjef&ji 


®j/0j 


Equations  52  and  53  are  the  invariant  conditions  for  the  jet-momentum 

coefficient.  Note  that  if  a is  selected  to  be  other  than  6.  or  0. 

J J 
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(with  a selected  correspondingly),  Eq.  53  permits  three  optional  choices 
for  fi/fi,  but  if  o is  selected  to  be  6^  or  0.  (with  a selected  corre- 
spondingly) only  two  options  are  available  for 


The  relation  between  the  pressure  coefficients  at  corresponding  points 
in  the  two  spaces  is  given  by 

- 4-  <4§>  - <-^->  <-  4-  ¥■■> 


p u '3x'  vsusx'  ' U<x>  9x' 


or 


(— ^)  Cr 


S„s 


u°x 


(54) 


The  factor  (s^/s^)  in  Eq.  54  can  be  written  as  (s^/susx)2/3(s^/susx)1/3 
and  appropriately  manipulated  to  give 


St,  a 2/3  e i/3 
(71“)  “ {§>  & 


(55) 


Substituting  Eq  55  in  Eq.  54  and  indicating  the  functional  dependence  of 
the  pressure  coefficient  yields 


(e/fi2)1/3cp  [x*.  y*.  (Efi)1/3z*;  IG,  Ip] 
= (e/fi2)1/3Cp  |i*,  y*,  (Efi)1/35*;  1G,  IpJ 


where 

xg  = y*i»  y*o»  xc»  xt.  XK  (57) 

Ip  H ((tan  A) /B J , BJV,  Tr/a,  Kr/o,  a/a,  (fl/a)w,  (efi) 1/3 (K^Cj  )(58) 


and  is  given  by  Eqs.  9 and  10.  The  parameter  a may  be  selected  as 
any  one  of  the  parameters  xr,  Kr,  a,  0 ^ , or  <5^  depending  upon  the  type  of 
problem  under  consideration.  The  form  of  the  parameter  e (and  hence  e) 
is  arbitrary  for  subsonic  linear  flow,  and,  in  accordance  with  Eq.  20, 
is  e = = GJ)  for  nonlinear  transonic  flow.  The  optional  choices 

for  the  form  of  n (and  hence  h)  are  specified  by  Eq.  53.  The  value  of  u 
is  arbitrary  and  depends  upon  the  type  of  problem  under  consideration. 
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Common  selections  for  oj  are  (-1)  and  (2/3)  respectively  for  subsonic 

and  transonic  flow.  When  a is  selected  to  be  0j  or  S^,  the  underlined 

parameter  in  Eq.  58  is  omitted.  Additionally,  a selection  for  a that 

yields  unity  for  any  one  of  the  ratioed  parameters  in  Eq.  58  indicates 

a lack  of  an  invariance  constraint  on  that  ratio.  The  invariants  I„  and 

b 

Ip  are  separately  identified  since  the  Ig  invariants  are  purely 
geometric  at  all  times,  whereas  the  Ip  quantities  may  consist  of  a mixture 
of  geometric  and  fluidynamic  invariants,  the  exact  nature  of  which 
depends  upon  the  specific  forms  selected  for  a,  ft,  and  e.  There  is  no 
unique  combination  of  the  similarity  parameters  appearing  in  Eqs.  56  and 
58.  The  parameters  shown  may  be  rearranged  to  give  a wide  variety  of 
combinations  by  appropriate  manipulation  of  Eqs.  34,  36,  37  and  53 
(noting  the  equivalence  of  the  ratios  on  the  right-hand  side  of  Eq.  53.) 
Such  a procedure  is  permissible  providing  the  total  number  of  invariant 
parameters  remains  unchanged  in  any  rearrangement. 

The  following  aerodynamic  coefficients  are  defined 
CL  = cLr  + cLr 

► 

<CP£  - Cpu)dx*dy* 


% 
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"V  =_  2S~  x*  (Cp  - CPu)dx*dy* 


« I 


1 

nR  = 2S~  X*j  Cj  0j  ^)dy* 


where  (-0^ ) is  an  appropriate  downwash  angle  and  the  small -angle 

assumption  (sin0.=0.)  has  been  employed  in  the  jet-reaction  coefficients. 
J J 

Application  of  the  similarity  rules  yields  the  following  results 

Cp(x*,  y*,  0)  = (|-)1/3P  (x*,  y*.  0;  IG,  If)  (66) 


cLr  = dG.  If) 

(67) 

cLr  = (^-)1/3LR  (!G*  !f) 

(68) 

Cl  = (fV/3l  (IG.  Ip) 

(69) 

cDl  = C~)  1/3Pi  de.  If) 

(70) 

n2  1/3 

Cmr  = <I“>  Mr  (IG.  If) 

(71) 

CnR  » dG.  *F> 


Cm  = ttG.  XF> 

where  Iq  and  Ip  are  given  by  Eqs.  57  and  58,  respectively,  and 


L = Lp  + LR  and  M = Mp  + MR. 
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Since  in  the  linear  case  the  form  of  c may  be  selected  arbitrarily, 

it  is  of  interest  to  examine  the  consequences  of  several  choices  for  this 

2 

parameter.  For  e = £1  (e  = £2  ),  Eq.56  shows  that  the  pressure  coefficients 

at  affinely  related  points  in  the  original  and  transformed  spaces  are 

3 

identical.  For  e = 6^  , Eq.  25  shows  that  £1=0,  from  which  it  follows 
from  Eq.  41  that  a = a and,  from  Eq.  45,  that  0.  = 0..  Finally, 

J J 

equations  34,  35,  and  36  show  that  the  thickness  ratios,  camber  ratios, 
and  angles  of  attack  are  identical  in  the  original  and  transformed  spaces. 
For  e = IU<,  Eqs.  19  and  20  show  that  the  similarity  parameters  will  be 
identical  for  both  linear  and  nonlinear  flows.  This  allows  theoretical 
solutions  of  the  linear  and  nonlinear  flows  to  be  plotted  on  the  same 
set  of  graphs  in  terms  of  the  same  set  of  parameters.  The  two  theories 
would,  of  course,  yield  different  curves  on  such  plots.  The  curves  for 
linear  theory  would  be  valid  for  purely  subsonic  flows,  but  may  or  may 
not  be  valid  for  transonic  flows.  The  curves  for  nonlinear  transonic 
flow  would  be  valid  not  only  for  transonic,  but  for  subsonic  and 
supersonic  small -perturbation  flows  as  well. 

When  the  foregoing  similarity  rules  are  applied  to  two-dimensional 
flows,  the  finite-span  dependent  terms  are  obviously  omitted. 


AFFDL-TR-76-86 


SECTION  III 

COMPARISON  WITH  PREVIOUS  ANALYSES 

For  the  unblown  wing  there  is  no  dependence  upon  c.  and  6.,  and  the 

J J 

foregoing  results  may  be  shown  to  reduce  to  those  of  various  investigators 
by  appropriate  selection  of  ft,  e,  and  a,  and,  perhaps,  some  additional 
manipulation  of  the  invariants  as  previously  noted.  The  results  of 
Spreiter  (Appendix  B,  Ref.  5),  for  example,  may  be  obtained  by  using 

Eqs.  28  and  29  to  replace  A tan  A with  6^  tan  A and  taking  w=l , 

3 2 °° 

ft=B  /£,  e=8  A,  and  a=i  for  linear  flow,  and  taking  cj=1/3,  ft=x  , 

oo  oo  p r* 

e=kU  , and  a=r  for  nonlinear  flow. 

00  p 

As  mentioned  previously,  there  are  no  other  similarity  analyses  (to 
the  best  of  this  writer's  knowledge)  for  nonlinear  jet-flapped-wing  flow, 
hence  comparisons  will  be  limited  to  linear  flow.  Comparisons  will  be 
made  with  the  work  of  Siestrunk  (Ref.  1),  Levinsky  (Ref.  2)  and  Elzweig 
(Ref.  3).  Each  of  these  investigators  uses  a different  scaling  law. 

There  are,  perhaps,  other  works,  but  the  foregoing  provide  sufficient 
variety  for  the  present  purposes.  Among  the  aforementioned  analyses, 
Levinsky's  is  the  only  one  applicable  to  finite-span  wings;  the  other 
treatments  are  for  two-dimensional  flow.  None  of  the  analyses  include 
second-order  terms  in  the  jet  internal  flow;  hence,  in  comparing  with 
the  present  work,  the  parameter  will  be  taken  as  unity. 

Siestrunk  employs  a scaling  that  maintains  equal  angles  of  attack 
and  jet  deflection  angles  in  the  original  and  transformed  spaces. 

Levinsky  employs  Gothert's  scaling  (see  Ref.  6),  which  is  applicable  also 
to  axially  symmetric  flow.  Elzweig  employs  a scaling  for  which  the 
wing-surface  pressure  coefficients  are  equal  at  affinely  related  points 
in  the  original  and  transformed  spaces. 

The  results  of  Siestrunk  are  obtained  by  taking  t =k  =o,  a=a,  e-6  3, 

i i 00 

ft=l,  and  u)=-l.  Those  of  Levinsky  are  obtained  by  taking  t =k  =o,  a=a, 

/I 

e=8  ,ft=8  - , u)=-l.  The  results  of  Elzweig  are  obtained  by  taking 

2 

t^k^o,  a= a,  e=$m  , 0*$^  and  of-1. 
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SECTION  IV 

COMPARISON  WITH  EXPERIMENT 

As  noted  previously,  for  jet-flapped  wings  there  is  no  experimental 
data  of  appropriate  configurational  simplicity  and  of  sufficient  gener- 
ality in  parametric  variations  to  permit  validation  of  the  nonlinear 
rules  (to  the  knowledge  of  this  writer).  The  data  available  for 
validation  of  the  linear  rules  is  rather-  sparse  and  is  plagued  with 
uncertainties  regarding  wind-tunnel  wall  corrections.  Nevertheless  some 
limited  comparisons  are  made  in  this  section.  The  data  sources  selected 
for  comparison  are:  (1)  Air  Force/Northrop  tests  (Refs.  8 and  9), 

(2)  Air  Force/Convair/Canadian  tests  (Refs.  10,  11,  12,  and  13),  and 

(3)  French  O.N.E.R.A  tests  (Ref.  14).  The  first  two  test  series  are 
for  two-dimensional  airfoils,  whereas  the  third  is  for  a finite-aspect- 
ratio  wing  employing  a semi  span  test  arrangement. 

AIR  FORCE/NORTHROP  TESTS 

These  tests  (Refs.  8 and  9)  were  conducted  at  the  Arnold  Engineering 
Development  Center  (AEDC)  in  the  4T  wind  tunnel  which  has  a 4-x  4-foot 
porous-wall  test  section.  The  model  tested  has  a modified  NACA  64A406 
airfoil  section,  a 10-inch  chord,  a 20-inch  span  and  was  mounted  between 
two  large  end  plates.  Tests  were  conducted  at  Mach  numbers  of  0.70,  0.80, 

0.85,  0.90  and  0.95  and  at  effective  jet-flap  deflection  angles  (S.)  of 

J 

0,  35,  55,  and  88  degrees.  Forces  and  moments  were  not  measured  directly, 
but  were  deduced  from  measured  surface  pressures.  The  data  presented 

4 

in  Refs.  8 and  9 was  not  corrected  for  tunnel  or  end-plate  interference 
effects.  Therefore,  for  the  present  comparisons  an  approximate 
correction  is  derived  in  Appendix  B of  this  report. 

Since  the  camber  line  of  the  airfoil  used  in  these  tests  is  almost 
parabolic,  the  parabolic  camber-line  relations  of  Appendix  C are  used  in 
estimating  the  airfoil  aerodynamics.  The  applicable  incompressible  lift 

I 

4 

Private  communication  with  the  test  conductor. 
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coefficient  relation  is  given  by  Eq.  5 in  Appendix  C.  Taking  e = g^3 * 5, 
fi  = 1 , a = a,  and  u>  = -1  in  Eq.  69  yields,  for  the  subsonic  similarity 
relation, 

[^K  (e“KlrCj)  + aC^a  (6»KltCJ)  + 63ct6l  W]  (74) 

The  consequences  of  Eq.  74  are  compared  with  experiment  in  Fig.  3. 

The  35-degree  jet-flap  deflection  case  was  selected  for  comparison  since 
the  theory  is  more  strictly  applicable  to  smaller  jet-flap  angles.  The 
experimental  data  shown  is  for  a constant  geometrical  angle  of  attack  , 
which  differs  from  the  actual  aerodynamic  angle  of  attack  due  to 
tunnel  interference  effects.  It  would  have  been  preferable  to  present 
the  experimental  results  at  a constant  aerodynamic  angle  of  attack,  but 
the  data  of  Refs.  8 and  9 for  = 35°  was  too  meager^  to  permit 
construction  of  the  required  cross  plots.  The  theoretical  results 
directly  comparable  to  experiment  are  shown  by  the  solid  symbols.  These 
results  were  calculated  using  the  aerodynamic  angles  of  attack  obtained  by 
the  method  of  Appendix  B.  The  corresponding  aerodynamic  angles  of 
attack  are  listed  in  the  table  on  the  figure.  Point  comparisons  are 
not  shown  at  supercritical  Mach  numbers  since  the  tunnel  interference 
analysis  is  not  valid  in  this  regime.  The  theoretical  lift-coefficient 
variation  with  Mach  number  at  a constant  angle  of  attack  is  shown  by  the 
solid  lines  in  Fig.  3.  The  constant  angle  of  attack  is  the  average  of 
the  aerodynamic  angles  of  attack  for  M^O.7  and  0.8.  These  limited 
comparisons  show  that  fairly  good  agreement  is  achieved  between  theory 
and  experiment. 

The  relative  importance  of  the  jet-sheet  second  order  effect  can  be 
seen  by  comparing  the  dashed-line  curve  for  K =1.0  with  the  solid  line 


3The  angle  of  attack  here  is  relative  to  a chord  line  through  the  leading 

and  trailing  edges  of  the  mean  camber  line.  The  angle  of  attack  (say  a) 
employed  in  Refs.  8 and  9 is  relative  to  an  arbitrary  chord  line.  The 
two  angles  of  attack  are  related  by  ot  = a - 0.60. 

5 

The  bulk  of  the  data  presented  was  for  6j  = 88°,  an  unacceptably  high 
value  with  which  to  make  valid  theoretical  comparisons. 
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curve  for  the  variable  shown.  In  this  case  the  second  order  effect 
is  significant.  The  dashed-dot  line  in  Fig.  3 shows  the  incremental  lift 
coefficient  due  to  camber  with  blowing,  (Ac£)k,  where 


! t [%«.Vj>-=<.,  <°>] 


Details  regarding  the  calculation  of  C£k  are  given  in  Appendix  C.  The 
incremental  bl own-camber- 1 ine  contribution  is  frequently  neglected  for  the 
reasons  noted  in  Appendix  C.  For  the  example  of  Fig.  3,  however, 
although  (Ac^  is  small,  it  is  of  sufficient  magnitude  to  warrant  its 
inclusion. 


AIR  FORCE/CONVAIR/CANADIAN  TESTS 

These  tests  (Refs.  10,  11,  12,  and  1j)  were  conducted  in  the 
Canadian  National  Aeronautical  Establishment  (NAE)  two-dimensional, 
high-Reynolds  number,  transonic  wind  tunnel  at  Ottawa,  Ontario.  This 
tunnel  has  a 60x60-inch  test  section,  porous  upper  and  lower  walls,  and 
employs  sidewall  inserts  to  form  a 15x60-inch,  two-dimensional  test 
section.  Several  different  jet-flapped  supercritical -type  airfoils  were 
tested  in  the  series  of  tests  reported  in  the  cited  references.  The 
airfoil  with  which  comparisons  are  made  herein  was  designated  as  an 
NAE  001002  airfoil,  possessed  considerable  aft  camber,  and  had  a chord 
of  15  inches.  Forces  and  moments  were  measured  by  sidewall  balances  with 
supplementary  data  taken  by  a wake  survey  rake.  A tunnel -interference 
angle-of-attack  correction  of  the  same  form  as  Eq.  1 in  Appendix  B was 
applied.  The  factor  k,  however,  was  obtained  from  drag  and  axial 
momentum  considerations  rather  than  in  the  manner  of  Appendix  B. 
Considerable  spread  in  the  value  of  k was  found,  but  a single  mean  value 
of  0.55  per  degree  was  employed  in  the  data  reduction. 

A comparison  of  the  results  given  by  the  subsonic  similarity  relation 
of  Eq.  74  with  some  NAE  test  data  is  given  in  Fig.  4.  Figure  4 shows 
that  rather  good  agreement  is  obtained  for  subcritical  Mach  numbers. 

Note  also  that  the  lift-coefficient  increment  due  to  the  blown  camber 
line  is  significant. 
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Some  comment  is  appropriate  regarding  the  reference  angle  of 
attack,  a1,  employed  in  Fig.  4.  When  the  unblown  experimental  data  used 
for  Fig.  4 was  reduced  to  zero  Mach  number  by  the  Prandtl-Glauert 
relation,  the  angle  of  zero  lift  was  found  to  occur  at  a 0.66-degree 
angle  higher  than  predicted  by  thin-airfoil  theory.  This  discrepancy 
is  assumed  to  be  due  to  a constant  geometric  error  in  the  experiments, 
which  is  equivalent  to  considering  the  angle  of  attack  to  be  measured 
from  a chord  line  different  from  that  specified.7  In  the  present  case, 
the  angles  of  attack  (in  degrees)  relative  to  the  specified  and  apparent 
chord  lines  are  related  by  a’  = a - 0.66. 

The  camber-line  aerodynamics  for  Fig.  4 was  calculated  by  the  method 
of  Appendix  C for  an  arbitrary  camber-line  shape.  Note  also  that  the 
angles  of  attack  (in  degrees)  relative  to  the  specified  chord  line  and  the 
camber-line  chord  are  related  by  a = a - 0.34. 

The  relative  importance  of  the  jet-sheet  second  order  effect  can  be 

seen  by  comparing  the  dashed-line  curve  for  = 1.0  with  the  solid-line 

curve  for  the  variable  K shown.  In  this  case  the  second-order  effect  is 

n 

negligible.  For  a jet-momentum  coefficient  of  a magnitude  sufficient 
to  choke  the  jet  nozzle  at  a free-stream  Mach  number  less  than  0.5,  the 
second-order  effect  would  be  significant. 

FRENCH  O.N.E.R.A.  TESTS 

The  French  tests  (Ref.  14)  were  conducted  on  a half-span  rectangular 
wing  model  of  aspect  ratio  3.4.  The  airfoil  section  is  an  NACA  64A010 
airfoil  truncated  and  modified  at  the  88%  chord  position  yielding  a 
11.4%  thick  airfoil.  Further  details  on  the  model  and  facility  were 
not  available. 


7This  assumption  is  justified  by  an  improved  agreement  between  theory 
and  experiment. 
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In  applying  the  subsonic  similarity  law,  Hartunian's  theory  (Ref.  15) 
for  a finite-span  jet-flapped  wing  is  employed  to  obtain  the  wing 
characteristics  at  zero  Mach  number.  On  the  basis  of  Hartunian's 
relations,  the  lift,  induced  angle  of  attack,  and  induced  drag, 
respectively,  are  given  by: 

“I.  fc.i>  * 

(Cl)Mo=0  (Cj)  (76) 

i + 5j 

irA  + 2cj 


(CL)Ma)=0 
(_0i)Mco=O  = irA  + 2c, 


( Di)M0O=0  -0  (CL)M0O=0  tiA  + 2cj 


Taking  £=6^  , fl-l , a=a,  and  w=-l  in  Eq.  69  yields  the  similarity  relations 


„ 1_  [ac^(^Vj>  + 5j  'S  (B-Vj)] 


C*«j 

+ 26a)Kwcj 


(-0±)  = (-0i)»  =0 


% = (~0i) 


M„=0  CL 


Determination  of  in  Eq.  79  requires  values  for  the  relative 

nozzle  height,  h./c,  and  the  nozzle-flow  coefficient,  c . Since  neither 
j n 

of  these  quantities  was  available  in  the  referenced  document,  guessed 
values  were  employed.  For  this  purpose,  the  nozzle-flow  coefficient 
was  taken  to  be  unity  and  hj/c  to  be  0.008.  The  value  of  h^/c  was 
selected  in  the  following  manner.  For  tests  conducted  on  an  NACA  0018 
airfoil  in  Ref.  14,  the  average  value  of  hj/c  was  0.0129  (see  Fig.  10 
of  Ref.  14).  Ratioing  this  value  according  to  wing  thickness  yields 

hj/c  = 0.008. 
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A comparison  of  the  lift-coefficient  scaling  given  by  Eq.  79  with 
the  experiments  of  Ref.  14  appears  in  Fig.  5,  where  reasonably  good 
agreement  is  obtained  for  a jet-deflection  angle  of  29°,  but  poorer 
agreement  is  obtained  for  6j=-2°.  Apparently  the  influence  of  is 
negligible  in  this  case. 
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SECTION  V 
CONCLUSIONS 

Applying  the  known  principles  of  scaling,  linear-subsonic  and 
nonlinear-transonic  similarity  rules  have  been  derived  for  a finite- 
span-jet-flapped  wing  with  partial  or  full-span  blowing.  In  deriving 
the  rules,  an  attempt  was  made  to  keep  the  presentation  as  general  as 
possible  in  order  to  display  the  interrelation  between  the  linear  and 
nonlinear  rules  and  to  allow  freedom  in  adjusting  the  scaling  to 
emphasize  the  parameters  and  type  of  scaling  most  appropriate  to  a 
particular  problem  under  consideration.  Although  this  generalized 
approach  makes  the  analysis  slightly  more  difficult  to  follow,  it  is 
believed  to  be  worthwhile  for  the  additional  visibility  it  provides. 

The  effect  of  jet-supply  pressure  ratio  was  delineated  by  considering 
second-order  quantities  in  the  jet-sheet  compatibility  condition.  The 
importance  of  jet-supply  pressure  ratio  depends  upon  the  flight  Mach 
number  at  which  the  jet  nozzle  (assuming  a convergent  configuration) 
chokes.  At  this  time  it  is  not  known  whether  typical  flight  vehicles 
will  fall  within  the  parametric  spectrum  where  jet-supply  pressure  ratio 
is  an  important  consideration.  In  comparisons  with  three  different  sets 
of  wind-tunnel  data,  jet-supply  pressure  ratio  was  found  to  be  of 
significant  importance  for  only  one  set  of  data  considered. 

This  investigation  also  disclosed  that  camber-line  effects  with 
blowing  assume  more  importance  than  heretofore  was  the  case  in  jet- 
augmented  mechanical  flap  applications. 

Finally,  although  reasonably  good  agreement  was  obtained  for 
limited  comparisons  with  experiment  for  linear-subsonic  flow,  there  is 
a critical  need  for  well-designed  wind-tunnel  experiments  to  validate 
both  the  linear  subsonic  and  nonlinear  transonic  similarity  rules.  The 
experiments  should  be  planned  to  eliminate  tunnel  interference  effects 
or  should  be  conducted  in  a facility  where  the  magnitude  of  the  inter- 
ference can  be  accurately  predicted.  In  the  tests  parametric  variations 
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should  be  made  within  ranges  consistent  with  the  small  perturbation 
assumption  of  the  similarity  analysis.  Parametric  variations  also 
should  be  specifically  tailored  to  test  three-dimensional  effects  and 
the  nonlinear  transonic  rules. 
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APPENDIX  A 

SECOND-ORDER  JET-SHEET  COMPATIBILITY  CONDITION 
FOR  COMPRESSIBLE  FLOW 


For  the  jet  sheet  to  be  compatible  with  the  external  flow,  the 
jet-sheet  internal  static  pressures  at  the  sheet  upper  and  lower 
boundaries  must  be  equal  to  the  corresponding  static  pressures  in  the 
external  flow,  and  the  sheet  boundaries  must  be  stream  surfaces  of  the 
external  flow.  In  deriving  the  compatibility  condition,  the  lateral- 
flow  velocity  perturbations  are  assumed  to  be  negligible  compared  to 
the  vertical  and  longitudinal  ones,  thereby  permitting  the  jet-sheet 
flow  to  be  treated  as  two-dimensional  in  any  plane  where  y = a constant. 

A similar  approach  is  employed  by  Maskell  and  Spence  in  their  treatment 
(Ref.  16)  of  a finite-span  jet-flapped  wing.  With  reference  to  Fig.  2, 
it  is  also  assumed  that  the  jet-center-line  radius  of  curvature,R,  is 
large,  the  jet  thickness,  h,  is  very  small,  and  the  downward  displacement 
of  the  jet  is  small,  such  that  h/R  « 1 and  R~^  = - 30j$/3x. 

In  jet-flap  applications  the  jet-exhaust-nozzle  height,  h.,  is 

J 

usually  rather  small  due  to  the  geometrical  constraints  imposed  by  the 
thinness  of  the  airfoil  trailing  edge.  Consequently,  for  a convergent 
nozzle,  choked  (critical)  flow  may  occur  if  the  jet-supply  pressure  is 
sufficiently  high.  Assuming  the  jet-sheet  internal  flow  is  isentropic1, 
the  jet  flow  downstream  of  the  nozzle  will  have  a differing  character 
for  subcritical  and  supercritical  nozzle  flow  conditions.  For  sub- 
critical  nozzle  flows  the  magnitude  of  the  nozzle  exit  pressure  will  be 
essentially  governed  by  the  external  freestream  static  pressure.  For 
the  limiting  case  of  zero  jet  thickness  in  Spence's  incompressible  flow 
analysis  (Ref.  17),  it  is  implicit  that  the  static  pressure  along  the 
jet  center  line  is  equal  to  the  freestream  static  pressure.  As  will  be 
seen  subsequently  in  this  appendix,  this  condition  is  also  approximately 


In  the  real  flow,  shock  waves  are  likely  to  be  present.  The  assumption 
of  isentropy  is  an  approximation,  strictly  valid  for  supercritical  duct- 
supply  pressures  only  slightly  exceeding  the  critical  presssure 
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true  for  a first-order  compressible  flow.  If,  however,  the  nozzle  flow 
is  supercritical,  the  nozzle  static  pressure  is  governed  by  the  jet-supply 
pressure  and  may  on  the  average  be  greater  than  the  freestream  static 
pressure.  In  the  analytical  modeling  for  this  case,  therefore,  some 
provisions  must  be  made  for  accounting  for  the  longitudinal  decay  of 
the  center-line  over-pressure  within  the  jet.  This  can  be  accomplished 
by  a second  (or  higher)  order  analysis  of  the  jet-sheet  internal  flow, 
as  in  this  appendix. 


With  reference  to  Fig.  2,  the  internal  and  external -flow  field 
horizontal-velocity  components  are  given  respectively  by 


u = U™,  + u' 


u = U + u' 


( A— 1 ) , (A-2) 


where  u1  « U and  u'  « U 


The  irrotationality  condition,  in  natural  coordinates  (n,s),  is  given  by 
Ref.  6 as 


1 3u  _ l 
~ 3n  r 


(A-3) 


Expanding  r along  n gives 


r = R + (3r/3n)  n + 
o 


(A-4) 


Substituting  Eq.  A-4  in  Eq.  A-3,  integrating,  and  applying  the  condition 


u=uQ  at  n=0,  yields 


u exp 
iL 


n _ 1 ,3r.  n2  , _ 


R 2 3n  0 r2 


(A-5) 


Assuming  [l-(3r/3n)Q]«  1 and  (n/R)  « 1,  making  use  of  Eq.  A-l,  and 
expanding  Eq.  A-5  yields  to  second  order  in  small  quantities 


u'  = S’o  + + u'o)  (n/R) 


(A-6) 
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P 


To  an  order  consistent  with  nonlinear  transonic-flow  theory,  the 
pressure  in  the  external  stream  is  given  by 

P = Pc  ~ Poouco2  (u’/uj  (A 

The  jet-sheet  internal  pressure  to  second  order  is 


P ~ Poo  “ p, 


Along  the  center  line,  Eq.  A-8  gives 


?o  = Poo  - p JJ  + \ d-H.2) 

U„  1 Uoo2 


The  internal  pressures  at  the  upper  and  lower  jet  boundaries  are  given 
by  substituting  u'y  and  u'£  respectively  in  Eq.  A-8,  where  from  Eq.  A-6, 
u‘u  and  u'^  are  given  respectively  by 

«'u  * “'o  +1  (Uoo  + u'oXh/R)  (A-10) 

u’n  « u’o  - I (u~  + u'oHh/m  (A-l  1 ) 


Applying  the  boundary  conditions 


Pu  ” Pu  and  Pe  “ P£ 


(A-l 2) 


and  employing  Eqs.  A-8,  A-10  and  A- IT  yields  to  second  order,  the  result 


Uoo2  <£>  [l  + (2-M,2)  (Jr*)]  (A-l 3) 


Pu  - Pe  = - Poo 


The  variable  h ( s ) may  be  determined  from  global  continuity.  Application 
of  continuity  between  downstream  infinity  and  an  arbitrary  s-station 
yields 


(p/p«,)  (u/UoJdn 


(A-l 4) 
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Employing  Eq.  39  of  Ref.  4,  the  density  to  second  order  is 

Making  use  of  Eqs.  A-2,  A-6,  and  A-15  in  Eq.  A-14  yields,  to  second 
order 

* - [i  - d-«.2)  hr  (#-’6) 

Uoo 


Substituting  Eq.  A- 1 6 in  Eq.  A-13  yields,  to  second  order 


Pu  " Pt  = -(P„A,hoo>  (!+  jp  ) r* 


(A— 1 7) 


Solving  Eq.  A-9  for  (u^/U^)  and  substituting  the  result  in  Eq.  A-17 
yields,  to  second  order 


Pi  “ Pu  — cu 


(A-18) 


where 


C.51- 


Po  - p„ 


’p  M 

<»  UO 


(A-19) 


and  R has*  been  replaced  by  the  approximation  R"^=  -80is/3x.  The 
parameter  p,  which  varies  longitudinally,  is  a factor  modifying  the 
usual  jet-sheet  compatibility  condition. 


It  is  of  interest  to  examine  C in  more  detail. 

7T 


If,  from  Eqs.  A-8,  A-10,  and  A-ll,  an  expression  for  (py  + pf)  is 
derived,  and  in  this  expression  Eqs.  A-9,  A-16,  and  A-12  are  substituted 

A 

respectively  for  (uQ/Uj,  (h/R),  p , and  p^,  the  following  result  is 
obtained  to  second  order 


Pu  + Pfc  , Poo  Ha 


(A-20) 


28 


AFFDL-TR-76-86 


In  the  vicinity  of  the  jet  sheet,  let 


+ u' 


(A-21 ) 


where  u'  and  u*  are  respectively  the  symmetrical  and  anti  symmetrical 

S a S 

parts  of  the  perturbation  velocities  relative  to  the  jet  center  such  that 


's»  = u * s * (u'as)u  = 


" <u’aS)i 


(A-22),  (A-23) 


From  Eqs.  A-7,  A-22,  and  A-23,  it  is  easily  shown  that 

(1/2)  (Pu  + P*)  = P.  - P^2  (u'g/OJ  (A-24) 

Recalling  that  p^  = p^,  Eqs.  A-20  and  A-24  yield,  to  first  order 


P0  ~ P» 

YPcoHx,2 


(A-25) 


The  velocity  u's  is  a free-stream  perturbation  velocity  at  the  jet 
boundaries  due  to  jet  thickness.  The  important  point  here  is  that  to 
first  order  PQ(s)  is  equivalent  to  the  pressure  distribution  along  the 
center  line  of  a thin  jet  of  finite  thickness  exhausting  at  zero 
inclination  into  a surrounding  stream.  Consequently,  pQ(x)  can  be 
determined  independently  of  the  more  general  curved  jet-sheet  flow 
problem. 


For  subcritical  nozzle  flow  into  a subsonic  external  stream  it  is 
well  known  that  p.=p  and  hence  it  is  reasonable  to  take  p (x)=p  for 
this  case.  For  choked  nozzle  flow  into  an  incompressible  external 
stream,  an  approxiate  solution  for  pQ(x)  could  possibly  be  found.  With 
pQ(x)  then  known,  Spence's  (Ref.  17)  integral  equations  could  be 
appropriately  modified  and  possibly  solved.  The  resulting  solutions 
could  be  extended  to  subsonic  linear  compressible  flow  by  means  of  the 
similarity  relations  herein.  The  solutions  so  obtained  would  be 

A 

functions  of  the  supply  duct  pressur-  ratio  (P/pm)  in  addition  to  c, 

J • 

Note,  however,  that  the  present  analysis  is  a small  perturbation  one, 
for  which  (p0-P00)/yp,1iM„2  « 1.0,  from  which  it  follows  that  -*■  1.0. 
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With  this  in  mind,  it  hardly  seems  worthwhile  to  treat  as  a variable. 
Instead,  a constant  average  value  <CiT>  will  be  sought  such  that 


<C  > t 1.0 

7T 


where 


CP/ ~ ) = (1/2  + y/2) 


+ y/2)y/y-1  = l.f 


<P0>  = 1/2  (Pj  + pj  = 1/2  (p*  + p ) 


(A-26) 

(A-27) 

= 1.4) 

(A-28) 

;,  assume  that 

t + pJ 

(A-29) 

This  yields 


«V  s 1 - 


<P*/Pj  - 1 

2>M  2 


(A-30) 


For  subcritical  nozzle  flow,  the  flow  field  across  the  nozzle  is 
affected  by  the  external  stream,  generally  is  nonuniform,  and  therefore 
is  inappropriate  for  use  as  a reference  stream.  In  the  previous  analyses, 
this  is  one  reason  for  referencing  the  local-flow  state  within  the  jet 
sheet  to  the  jet  internal  flow  at  downstream  infinity.  For  choked  or 
supersonic  nozzle  flow,  the  details  of  the  flow  across  the  nozzle  are 
governed  by  the  upstream  conditions  in  the  nozzle  duct,  and  the  duct 
contours  can  be  designed  to  give  a uniform  flow  at  the  nozzle  exit. 

Hence,  from  the  standpoint  of  flow  uniformity  the  nozzle  exit  flow  can  be 
used  as  a reference  stream.  However,  for  this  case,  when  the  nozzle  flow 
is  sonic,  the  assumption  of  small  perturbations  is  locally  invalid.  This, 
then,  is  a second  reason  for  employing  the  downstream  jet  flow  as  a 
reference  stream  in  the  previous  perturbation  analyses.  Tentatively, 
disregarding  the  second  reason,  a perturbation  analysis  from  a choked 
nozzle  stream  can  be  employed  to  yield  additional  insight  regarding 


A perturbation  analysis  paralleling  that  performed  for  will  now 
be  outlined  for  flow  perturbations  referenced  to  the  nozzle  exit. 
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For  convenience  in  discussion,  the  nozzle-exit  Mach  number  is  tentatively 
assumed  to  be  arbitrary.  For  the  aforementioned  analysis,  the  applicable 
equations  are  of  the  same  form  as  the  Eqs.  A-l  through  A-19,  with  the 
flow-state  variables  at  downstream  infinity  within  the  jet  sheet  replaced 
by  those  at  the  jet-nozzle  exit.  With  these  considerations  in  mind, 
evoking  continuity  within  the  jet,  and  taking  the  jet  exit  Mach  number 
to  be  unity,  the  following  result  may  be  deduced  by  analogy  with 
Eqs.  A- 1 8 and  A-19. 


P*-P„ 


30, 


p*  - pu  = - (p~u»hJuj,r4)a+  — -)  (^) 


YP. 


(A-31) 


Making  use  of  Eq.  38  in  Ref.  4,  it  may  be  shown  that  to  first  order 


a 

Hr  )=  1 


U 


oo 


P*_Pco 


(A-32) 


for  (P*-Pj/^P*  <<  1.0  substituting  Eq.  A-32  in  Eq.  A-31  yields  to 
second  order  the  result 


30. 


p.-p  = - C'  (pUh)U  (t— 

JC  ru  TT  oo  oo  oo7  on  'CW  7 


°°  3x 


(A-33) 


where 


(A-34) 


If,  as  before,  a mean  value  of  pQ  given  by  Eq.  A-29  is  employed,  the 


mean  value  of  C'  is 

IT 


For 


<C'  > = 1 - 

TT 


P*  - P« 
2yp* 


(A-35) 


P/P*  < (P/P.)  ^ 5 


(A-36) 


it  is  found  that 


C' 
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where  <CT[>  and  C'K>  differ  from  each  other  by  less  than  3%.  In  view  of 
its  simpler  form,  <C'it>  will  be  employed.  The  second  order  jet-sheet 


compatibility  condition  then  becomes 


-K^c  (-dr) 


K = 1.0 

Tt  » ,* 

1 P/P* 


(P/Poo  * P/P*) 
((P/p*)  < (P/p J < 5 
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APPENDIX  B 

INTERFERENCE  CORRECTION  FOR  AIR  FORCE/NORTHROP  TESTS 

The  small  magnitude  of  the  lift-curve  slopes  for  the  unblown-airfoil 
data  in  Refs.  8 and  9 at  the  sub-critical  Mach  numbers  of  0.7  and  0.8 
indicates  that  interference  effects  were  likely  present  in  the  tests. 
Although  there  is  no  independent  data  available  to  directly  support  this 
contention,  an  engineering  estimate  can  be  made  by  transforming  the 
slopes  to  the  incompressible-flow  regime  (Mm-0)  by  the  Prandtl-Glauert 
rule^  and  comparing  these  with  a variety  of  data  for  approximately 
similar  airfoils. 

Application  of  the  Prandtl-Glauert  rule  to  the  unblown  airfoil  data 
of  Figs.  A-l  and  A-2  in  Ref.  9 yields  incompressible-flow  lift-curve 
slopes  of  3.7  and  3.5  per  radian  respectively.  For  six  percent  thick, 
airfoils  of  the  NACA  6-series  type,  the  low  Mach  number  experimental 
lift-curve  slopes  documented  in  Ref.  19  for  eight  different  airfoils  are 
6.42,  6.30,  6.25,  6.25,  6.02,  6.02,  5.73,  and  6.19  per  radian.  The 
foregoing  comparisons  support  the  contention  that  the  data  of  Refs.  8 and 
9 a^e  in  error  due  to  interference  effects. 

The  average  lift-curve  slope  for  the  eight  airfoils  cited  above  is 
6.15  per  radian.  For  convenience  and  in  view  of  the  uncertainties 
involved,  however,  the  theoretical  value  of  2n  will  be  employed  in  the 
analysis  which  follows. 


‘ 


The  interference-induced  angle  of  attack  (geometrical  angle  of 
attack,  cig,  minus  the  actual  aerodynamic  angle  of  attack,  a)  in  the 
wind  tunnel  is  assumed  to  be  given  by 


1 


ag  - a = C£r/k 


(B-l) 


i 


i 


It  is  generally  conceded  (see,  e.g..  Ref.  18)  that  the  Prandtl-Glauert 
correction  to  the  lift-curve  slope  gives  reasonably  good  results  for 
subcritical  Mach  numbers. 
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where  k is  the  interference  factor  which  is  to  be  determined.  A similar 
correction  is  applied  to  the  data  of  Refs.  10  through  13,  except  that  k is 
determined  in  a different  manner.  Differentiating  and  manipulating 
Eq.  B-l  yields 

dcJ,r/da 

k = j-  8 

! _ dc*r/dag  (B-2) 

dcnr/da 


where  dc^^da^  is  the  lift  curve  slope  in  the  wind  tunnel  and  dcfc  /da 
is  the  true  lift  curve  slope  in  free  air.  Equation  B-2  applies  to  both  the 
unblown  and  blown  airfoils,  since  c£  is  the  circulation  lift  coefficient. 
However,  the  unblown  airfoil  data  only  is  used  to  obtain  k values  which 
are  then  applied  to  the  blown  airfoil  in  correlating  the  similarity  rules 
with  experiment.  The  lift-curve  slope  required  in  Eq.  B-2  is  obtained 
from 


dctr 

( -) 

v da  ;M 


. 3ctr 

6. 


3a  M =0 


where,  as  previously  stated,  the  incompressible-flow  slope  is  taken  to 
be  2ir. 


From  the  unblown  airfoil  data  of  Figs.  A-l  and  A-2  in  Ref.  9,  k was 
determined  to  be  12.38  and  13.41  per  radian  respectively  for  free-stream 
Mach  numbers  of  0.7  and  0.8. 
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APPENDIX  C 

ANALYSIS  OF  CAMBERED  JET-FLAPPED  AIRFOILS 

The  most  common  application  of  the  jet-flap  principle  to  date  has 
been  to  jet-augmented  mechanical  flaps  for  use  during  the  take-off  and 
landing  flight  phases.  For  this  application,  the  lift  coefficients 
achieved  are  very  high  and  the  relative  lift  contribution  due  to  camber- 
line jet-sheet  interaction  is  negligible.  For  pure  jet-flapped  wings  at 
high  speeds,  the  camber-line  contribution,  although  small,  can  be  a 
relatively  higher  fraction  of  the  total  lift  and  should  be  taken  into 
account  in  some  instances.  Because  of  previous  emphasis  on  very  high- 
lift  applications,  little  attention  has  been  given  to  the  analysis  of 
camber-line  effects.  The  sole  investigation  on  the  subject  appears  to  be 
that  of  Hough  (Ref.  20),  who  has  formulated  an  analysis  for  a polynomial 
camber  line,  but  has  provided  specific  numerical  results  only  for  the 
parabolic  case.  There  is  a need  for  a prediction  capability  for 
arbitrary  camber-line  shapes,  since  the  parabolic  camber  line  is  not 
necessarily  the  best  one  for  high-speed  applications.  Supercritical 
airfoils,  for  example,  are  possible  candidates  for  jet-flap  application. 
Many  of  these,  such  as  the  NAE  001002  airfoil  treated  herein,  are 
typified  by  considerable  aft  camber. 

Before  pursuing  the  subject  of  arbitrary  camber  further,  it  is  of 
interest  to  comment  on  Hough's  parabolic  results.  Hough  notes  that,  as 
a consequence  of  his  employment  of  a 3-control -point  calculation 
procedure,  his  results  may  be  in  error  at  low  values  of  the  momentum 
coefficient.  Since  at  high  speeds  the  momentum  coefficient  is  typically 
small  and  since  one  of  the  airfoils  examined  in  this  report  has  an  approxi- 
mately parabolic  camber  line,  it  was  deemed  of  interest  to  investigate 
the  magnitude  of  the  error  noted  by  Hough.  Consequently,  the  derivatives 
3c^/3k  and  3cm/3<  were  recalculated  using  a 9-control -point  calculation. 

The  resulting  curves,  along  with  a 3-point  calculation,  are  plotted  on 
Fig.  6.  Figure  6 shows  that  the  error  for  the  3-control -point  calculation 
increases  significantly  as  the  momentum  coefficient  approaches  zero. 
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One  could  further  develop  the  Hough  formulation  for  polynomial 
camber  lines  for  possible  application  to  arbitrary  camber  lines.  It 
appears,  however,  that  the  development  would  tend  to  become  unwieldly 
for  higher-degree  polynomials,  and  even  these  might  not  adequately 
represent  some  camber  lines.  Instead,  an  alternative  method  is  presented 
herein.  The  calculation  procedure  for  this  method  involves  only  a 
simple  quadrature,  is  easily  understood,  and  is  readily  performed  by 
hand. 

Applying  the  principle  of  superposition,  the  solution  for  a jet- 
augmented,  single-element,  mechanical-flapped  airfoil  given  in  Ref.  21 
can  be  used  to  obtain  a solution  for  the  segmented  airfoil  shown  in 
Fig.  7.  The  segmented  airfoil,  in  turn,  can  be  considered  to  be  an 
approximation  to  a continuous  camber  line.  On  this  basis,  the  lift 
coefficient  and  leading-edge  pitching  moment  coefficient  due  to  camber 
for  the  segmented  airfoil  of  Fig.  7 with  6^=0  are  given  respectively^  by 


(clK  = “ l <cj’  Xn)  (0n  * en-i> 

n=0 


= ~ I cm<s  (cjy  xn)  (On  ~ ®n-l^ 


n=0 


where  0, 


n 


yn+l  - yn 
xn+l  “ xn 


and  Xq  = 0 and  8_i 


(C-l) 


(C-2) 


(C-3) 


(C-4) 


A method  which  avoids  the  use  of  camber-line  slopes  and  is  more  accurate 
than  Eqs.  C-l  and  C-2  can  be  obtained  by  formulating  the  integrals 
resulting  from  taking  the  limiting  process  in  Eqs.  C-l  and  C-2.  These 
integrals  can  be  integrated  twice  by  parts  yielding  integrals  whose 
integrands  involve  the  camber-line  ordinate  multiplied  by  an  influence 
function.  The  influence  functions  required  in  this  approach  were  not 
available  at  the  time  of  this  investigation,  hence  Eqs.  C-l  and  C-2  are 
employed  in  the  present  calculations  as  a matter  of  expediency.  The 
details  and  influence  functions  for  this  second  method  will  be 
published  in  a technical  paper  at  a later  date. 
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The  derivatives  cAg  and  cm<5  are  given  in  Ref.  21  . In  applying  Eqs.  C-l 
through  C-4,  the  spacings  employed  can  be  adjusted  locally  according  to 
whether  0 is  varying  rapidly  or  slowly. 


A comparison  of  the  results,  using  Eq.  C-l,  with  a 9-control -point 
calculation  for  a parabolic  camber  line  is  shown  in  Fig.  6.  For  the 
calculation  shown  in  this  figure,  the  segmented  mean  line  consisted  of 
straight  lines  drawn  between  x stations  of  0,  0.05,  0.1,  0.2,  0.3,  0.4, 
0.5,  0.6,  0.7,  0.8,  0.9,  0.95  and  1.00  on  the  actual  mean  line.  As  can 
be  seen,  the  agreement  with  the  more  exact  calculation  is  excellent. 


In  general,  for  the  incompressible  flow  about  a pure  jet-flapped 
thin  airfoil,  the  lift  and  moment  coefficients  are  given  respectively 


where  the  derivatives  with  respect  to  a and  6.  are  given  in  Ref.  21  and 

J 

the  derivatives  with  respect  to  k for  a parabolic  camber  line  are  given 
in  Fig.  6. 


Additional  values  were  calculated  at  x=  .1,  .2,  .3,  .4, 

linear  interpolation  method  assuming  the  variation  with  x 
similar  to  the  c^=0  case. 
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Figure  3. 


Comparison  of  Theory  and  Experiment  for  Air  Force/NorthroD 

and^l”"  a-Mnd«t’nd  NACA  64A4°6  Airfoil  Section  (Refs.  8 P 
9J,  0.63  Deg,  6^  = 35  Deg,  6^/0=0.0020 
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Figure  4.  Comparison  of  Theory  and  Experiment  for  Air  Force/Convair/ 
Canadian  Tests  on  an  NAE  001002  Airfoil  Section  (Refs.  10, 
11,  12,  and  13);  a'  = -0.66  Deg,  = 30  Deg,  h./c=0.0030 
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Figure  6.  Comparison  of  Three  Different  Calculative  Procedures  for 
a Parabol ical ly  Cambered  Jet-Flapped  Airfoil 
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